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Navy  Case  No.  76257 

SYSTEM  AND  MEHTOD  FOR  RAPIDLY  TRACKING 
HIGHLY  DYNAMIC  VEHICLES 


j  STATEMENT  OF  GOVERNMENT  INTEREST 

i 

I  The  invention  described  herein  may  be  manufactured  by  or  for 

the  Government  of  the  United  States  of  America  for  Governmental 

li purposes  without  the  payment  of  any  royalties  thereon  or 

I 

I  therefor. 

I 

CROSS-REFERENCE  TO  RELATED  APPLICATIONS 
This  patent  application  is  co-pending  with  related  patent 
application  entitled  SYSTEM  AND  METHOD  FOR  RAPIDLY  TRACKING 
VEHICLES  OF  SPECIAL  UTILITY  IN  LOW  SIGNAL-TO-NOISE  ENVIRONMENTS, 
Navy  Case  No.  76256  by  the  same  inventors  as  this  application. 

BACKGROUND  OF  THE  INVENTION 
(1)  Field  of  the  Invention 

I  The  invention  relates  generally  to  the  field  of  signal 

II 

i  processing  and  more  particularly  to  a  system  and  method  for 
rapidly  detecting  a  moving  target  and  determining  its  movement 
characteristics,  such  as  range,  bearing,  speed  and  course  in  a 
noisy  environment. 


1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


(2)  Description  of  the  Prior  Art 

Detection  of  a  moving  object,  such  as  a  target,  and 
determination  of  its  range,  bearing,  speed  and  course  in  an  ocean 
environent,  is  a  difficult  task,  particularly  if  the  target  is 
moving  relatively  noiselessly  and  it  is  desired  to  perform  the 

detect  as  early  as  possible.  Typically,  acoustic  sensors  are 

j 

||  used  to  detect  acoustic  energy  (sound  waves)  emitted  by  a  moving 

j 

i object  and  convert  such  energy  to  electrical  signals,  and  complex 

j  j 

|l  signal  processing  operations  are  performed  in  connection  with  the 
electrical  signals  to  isolate  and  provide  the  desired 
information.  An  ocean  environment  is  generally  very  noisy,  and 
so  low-level  acoustic  signals  typical  of  quietly-moving  targets 
and  the  high  level  of  ambient  noise  joint  to  provide  a  relatively 
low  ratio  of  desired  signal-to-noise  in  the  electrical  signal 
provided  by  the  sensor,  which  makes  early  and  accurate  detection 
quite  difficult.  In  current  systems,  signals  that  do  not  have  a 
signal-to-noise  ratios  above  a  selected  predetermined  threshold 
value  are  ignored,  in  which  case  such  signals  are  not  available 

I  to  provide  information  which  may  potentially  be  useful  in 

i 

;  characterising  the  motion  of  the  target. 

t  i 

!  i 

I 

SUMMARY  OF  THE  INVENTION 

It  is  therefore  an  object  of  the  invention  to  provide  a  new 
and  improved  system  and  method  for  rapidly  tracking  moving 
objects  in  a  noisy  environment. 


2 


i  1 

In  brief  summary,  the  invention  provides  a  trajectory 

2 

estimation  system  for  estimating  a  trajectory  of  a  target  in 

3 

response  to  a  series  of  data  items  which  generated  in  response  to 

4 

motion  of  the  target.  The  trajectory  estimation  system  includes 

5 

a  data  segmentation  means  and  a  trajectory  selection  means.  The 

6 

data  segmentation  means  processes  the  series  of  data  items  in 

7 

accordance  with  a  regression/multiple-hypothesis  methodology  to 

8 

generate  a  plurality  of  segments,  each  having  associated  data 

9 

items  which  have  similar  features.  The  trajectory  selection 

10 

means  for  processing  said  segments  in  accordance  with  a  multiple- 

11 

model/hypothesis  methodology  to  generate  a  corresponding 

12 

statistically-supportable  candidate  trajectory  motion  estimate  of 

13 

i 

target  motion  thereby  to  provide  indicia  of  an  overall  trajectory 

14 

of  the  target. 

1 

15 

16 

BRIEF  DESCRIPTION  OF  THE  DRAWINGS 

17 

This  invention  is  pointed  out  with  particularity  in  the 

18 

appended  claims.  The  above  and  further  advantages  of  this 

19 

i 

invention  may  be  better  understood  by  referring  to  the  following 

20 

1  description  taken  in  conjunction  with  the  accompanying  drawings, 

21  i 

1  in  which: 

22 

FIGS.  lA  and  IB  together  constitute  a  functional  block 

23 

diagram  of  a  system  constructed  in  accordance  with  the  invention; 

24 

FIGS.  2A  through  3  comprise  flow  diagrams  illustrating  the 

25 

operation  of  the  system  depicted  in  FIGS.  lA  and  IB. 
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DESCRIPTION  OF  THE  PREFERRED  EMBODIMENT 
FIGS.  lA  and  IB  together  constitute  a  functional  block 
diagram  of  a  system  10  for  rapidly  tracking  highly  dynamic 
vehicles,  constructed  in  accordance  with  the  invention.  With 
reference  to  FIG.  lA,  the  system  10  includes  a  sensor  arrangement 
11  that  receives  acoustic  energy  (sound)  in  the  form  of  signals 

i  from,  for  example,  an  ocean  environment,  converts  the  signals  to 

j 

i electrical  form,  and  records  the  electrical  signals  for  later 

:  j 

I  processing.  A  fast  Fourier  transform  arrangement  12  performs  a 
conventional  fast  Fourier  transform  (FFT)  operation  in  connection 
with  the  recorded  signals  to  thereby  generate  phase  and  amplitude 
spectral  beam  maps  for  the  signals.  A  signal  follower  module  13 
receives  the  beam  maps  from  the  fast  Fourier  transform 
arrangement  12  for  signals  at  successive  points  in  time  and 
determines  whether  the  beam  map  indicates  that  the  signal-to- 
noise  ratio  of  the  signal  as  provided  by  the  sensors  11  exceeds  a 
predetermined  detection  threshold  value,  thereby  to  determine 
when  the  signals  represent  signals  from  a  particular  target  and 

i 

i  effectively  distinguising  such  target  signals  from  environmental 

I 

:  and  other  noise. 

i!  When  the  signal  follower  13  determines  that  a  beam  map  from 

I 

!  the  fast  Fourier  transform  arrrangement  12  exceeds  the 

i 

predetermined  detection  threshold  value,  a  detection  threshold 
comparator  14  compares  the  beam  map  corresponding  to  the  signal 
at  detection  to  the  beam  map  immediately  prior  to  detection  (that 
is,  for  the  last  beam  map  from  the  fast  Fourier  transform 
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arrangement  12  that  did  not  exceed  the  predetermined  detection 
threshold  value)  to  detect  similarities.  A  beam  regions  bound 
module  15  receives  the  beam  maps  and  similarity  information,  and 
bounds  the  beam  maps  based  on  a  priori  information,  such  as 
kinematic  and  other  information  known  about  likely  targets. 

The  detection  threshold  comparator  14  and  beam  regions  bound 

I  module  15  repeat  the  operations  with  each  beam  map  recorded  by 

I  j 

I I  the  sensor  arrangement  11  prior  to  the  signal  follower  13 

j  i 

ii  determining  that  a  signal  exceeded  the  signal-to-noise  threshold 
value.  This  allows  the  detection  threshold  comparator  14  and  the 
beam  regions  bound  module  15  to  obtain  information  concerning  the 
target  from  the  signals  recorded  prior  to  detection  (that  is, 
prior  to  the  signal  follower  module  13  determining  that  a  signal 
exceeded  the  signal-to-noise  threshold  value) ,  so  that  the  system 
10  will  not  have  to  rely  solely  on  signals  received  after  such 
time.  In  addition,  the  system  10  facilitates  a  restriction  on 
the  number  of  signals  that  it  will  have  to  analyze  and  allow  for 

j  subsequent  information  to  be  recorded  at  signal-to-noise  ratios 

I 

lower  than  the  detection  threshold  values.  In  particular,  a 

I 

I  measurement  track  formations  module  16  receives  the  information 

j 

from  the  beam  map  bounds  module  and  applies  a  lower  signal-to- 
noise  ratio  threshold  value  than  that  applied  by  the  signal 
follower  module  13  to  the  beam  maps  recorded  by  the  sensor 
arrangement  prior  to  signal  detection  as  determined  by  the  signal 
follower  module  13.  The  measurement  track  formations  module  16 
repeats  these  operations  through  a  series  of  iterations,  in  each 
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iteration  applying  a  lower  signal— to— noise  ratio  than  in  the 
previous  iteration,  to  extract  signal  information  from  the 
background  noise  and  clutter  in  those  beam  maps.  For  each  of  the 
beam  maps  that  satisfy  the  signal-to-noise  criteria  for  each  of 
the  iterations,  the  measurement  track  formations  module  16 
performs  an  inverse  fast  Fourier  transform  operation  to  transform 

the  bounded  beam  maps  to  provide  a  time-based  signal  for  later 

I 

I  processing. 

i 

I  The  signal  information  from  the  measurement  track  formations 

i 

module  16  is  then  used  by  a  data  segmentation  module  20  (FIG.  IB) 
and  a  trajectory  estimation  module  22  (FIG.  IB)  to  determine  the 
range,  bearing,  speed  and  course  of  the  target  which  is  the 
source  of  the  signal.  The  operations  of  the  data  segmentation 
module  20  and  the  trajectory  estimation  module  will  be  described 
below  in  detail  in  connection  with  FIGS.  2A  and  2B  (the  data 
segmentation  module  20)  and  FIG.  3  (the  trajectory  estimation 
module  22) .  Briefly,  however,  the  data  segmentation  module  20 
receives  the  signal  information  from  the  measurement  track 
formations  module  16  and,  using  that  information  and  a  priori 
kinematic  and  other  knowledge  concerning  likely  targets  from  a 

I  priori  knowledge  input  21,  generates  one  or  more  hypotheses 
regarding  movement  of  the  target.  The  trajectory  estimation 
module  22  receives  the  hypotheses  and  selects  one  as  the  most 
likely  hypothesis,  effectively  selecting  the  most  likely 
trajectory  (range  and  bearing)  of  the  target.  The  trajectory 
that  is  selected  is  verified  by  a  trajectory  analysis  and 
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validation  module  23  and  a  trajectory  characteristics  module  24 
using  conventional  statistical  measures  testing  the  likelihood  or 
probability  that  a  trajectory  is  representative  of  the 
information  contained  in  the  signals  received  by  the  sensor 
arrangement  11. 

As  noted  above,  the  data  segmentation  module  20  (FIG.  IB) 
generates  a  set  of  hypotheses  Hy  each  containing  one  or  more 

i 

I  segments  Sj.  Each  segment  Sj  is  a  hypothesized  line  segment  that 
the  data  segmentation  module  20  generates  in  response  to  the 
signal  information,  represented  by  a  series  of  data  items,  that 
the  data  segmentation  module  20  receives  from  the  measurement 
track  formations  module  16.  The  data  segmentation  module  20 
generates  the  segments  Sj  in  a  series  of  iterations  for  each 
successive  data  item  it  receives.  In  each  iteration,  the  data 
segmentation  module  20  effectively  attempts  to  add  the  data  item 
to  each  segment  3=  that  it  had  initiated  during  previous 
iterations,  and  generates  a  likelihood  measure  indicating  the 
likelihood  that  the  data  item  actually  belongs  to  each  of  the 
segments  Sj.  In  addition  during  each  iteration,  the  data 

!  segmentation  module  20  initiates  a  new  segment  3^  containing  only 
the  new  data  item,  for  the  possibility  that  the  data  item  is  the 
first  data  item  of  a  segment,  and  generates  a  likelihood  measure 
indicating  the  likelihood  that  the  data  item  is  the  first  data 
item  for  a  new  segment;  in  each  subsequent  iteration,  the  new 
segment  will  be  used  along  with  other  segments  initiated  during 
previous  iterations  as  possible  segments  for  subsequent  data 
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items.  In  addition,  during  each  iteration  the  data  segmentation 
module  generates  a  "false  alarm"  hypothesis  HjpA  for  the 
possibility  that  the  data  item  does  not  belong  to  any  segment. 

The  trajectory  estimation  module  22  prunes  the  hypotheses  and  the 
various  segments,  and  over  a  series  of  iterations,  the  data 
segmentation  module  20  and  trajectory  estimation  module  22 
I  cooperate  to  narrow  the  hypothesized  segments  Sj. 

The  operations  performed  by  the  data  segmentation  module  20 
land  the  trajectory  estimation  module  22,  each  during  one 
iteration,  are  depicted  in  FIGS.  2A,  2B  (data  segmentation  module 
20,  and  FIG.  3  (trajectory  estimation  module  22) .  With  reference 
initially  to  the  data  segmentation  module  20,  the  data 
segmentation  module  20  represents  each  segment  as  a  reduced  set 
of  regression  coefficients,  or  "features"  in  the  signal 
represented  by  the  data  stream.  With  further  reference  initially 
to  FIG.  2A. ,  upon  receiving  a  new  data  item,  identified  herein  as 
"Zj(t)"  (step  100),  the  data  segmentation  module  initially 
performs  a  series  of  steps  101  through  103  to  test  the 
I  statistical  consistency  of  the  data  item  Zi(t)  with  each  segment 

jiSj. 

ii  In  determining  the  statistical  consistency,  if  it  is  assumed 

j  that  a  segment  Sj  consists  of  "n"  data  items  previously  assigned 
I  to  the  segment  Sj,  the  data  segmentation  module  20  initially 
forms  a  predicted  data  item  value  Zi(t/n)  as  set  forth  in 
Equation  1; 


8 


1 


2^{t)  =  a^it*)  +  +V2a2(t*) 


( t-  t*) 


(1) 


2 

3 

4 

5 

6 
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ii  where  ai(t*)  are  the  i'th  time  derivatives  of  the  measurement 
I  Zi(t)  evaluated  at  time  t*  (step  101).  In  this  case,  the  data 
i|  segmentation  module  20  uses  t*  as  the  midpoint  of  the  successive 
i|  data  intervals  for  the  successive  data  items  to  minimize 

li 

!  I 

estimation  errors  for  the  regression  operation. 

After  generating  the  predicted  data  item  value  Zi(t/n)  via 
equation  (1) ,  the  data  segmentation  module  20  generates  a 
normalized  squared  residual  value  (step  102)  as 

f{t/n)  =  [Zi(t)  -  c{t/n)  [z^(t)  -2^{C/n)]  (2) 
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and  performs  a  chi-squared  test  in  connection  with  r(t/n)  to 
determine  whether  it  satifies  a  threshold  "gating"  value  (step 
I  103) ,  which  in  one  embodiment  is  set  to  the  value  thirty-six.  If 

i 

j  the  data  segmentation  module  20  makes  a  negative  determination  in 
ij  step  103,  it  will  ignore  the  data  item  received  in  step  100,  and 
will  returns  to  step  100  to  repeat  the  operations  in  steps  101 
through  103  in  connection  with  the  next  data  item. 

If,  on  the  other  hand,  the  data  segmentation  module  20  makes  a 
positive  determination  in  step  103,  it  assigns  a  probability 
value  Pa(Sj)  identifying  the  likelihood  that  the  new  data  item 


9 


belongs  to  segment  Sj  (step  104),  and  generates  an  updated 
segment  S/  including  the  new  data  item  (step  105)  .  The  data 
segmentation  module  2  0  generates  the  probability  value  Pa(Sj)  as  a 
function  of  the  normalized  squared  residual  value  generated  i 
step  102  prior  to  the  chi-squared  test,  in  a  conventional  manner. 
In  one  embodiment,  for  simplicity  the  probability  assignment  is 
I  obtained  by  mapping  the  normalized  squared  residual  (equation  2) 
i  in  an  intuitive  formula  that  approximates  the  complement  of  the 

i  I 

||  chi-squared  distribution  as  follows 


i 


1  -  ^  ^  for  f(C/n)  <  36 

36 

0  otherwise 


(3) 


The  data  segmentation  module  2  0  generates  the  updated  segment  Sj' 
to  include  the  new  data  item  (identified  as  the  '’n+l'‘st  data 
item)  in  the  following  manner.  Given  a  value  a(t„*)  as  the 
current  endpoint  of  segment  Sj,  the  data  segmentation  module  2  0 
generates  a  new  endpoint  a(t„+,*/n+l)  for  the  data  item  Zi(t)  as 

a{t^^^*/n+l)  =  a{t^^j^*/n)  +  K(n)  [z^(t)  -S^{t/n)]  (4) 


where 

a(t„^i*/i2)  =  t„*)  a{t^*/n)  (5) 


10 


and 


J 


^  ^n+X  ’ 


1 

0 

0 


0  1 


(6)  . 


The  new  endpoint,  together  with  the  points  previously  assigned  to 
the  segment  Sj,  defines  an  updated  segment  S/  .  The  data 
segmentation  module  20  determines  the  Kalman  gain  K(n)  for 
equation  3  as 

K{n)  ^  c{t/n)  (7) 

with 

c{C/n)  =  [H{C,  t„*)  R(t„*)  H{t,  +  0^]-^  (8) 

and 

H(t,  t^*)  =  [1  iC-t*)  ]  (9) 

and  the  corresponding  covariance  matrix  of  a(t„+i*/n+l)  is 

R{t„.^*/n+l)  =  [ J  -  K{n)A{t„,^*,  t„*)^]R{t„,^*/n)  (10) 

where 
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and  "I"  is  the  three-by-three  identity  matrix. 

After  generating  the  updated  segments  Sj'  for  all  segments 
for  which  the  chi-squared  test  was  satisfied  in  step  103,  the 
data  segmentation  module  20  effectively  updates  the  set  of 

( 

I  hypotheses  H;.  In  that  process,  the  data  segmentation  module  20 

ii  updates  hypotheses  Hjj  developed  during  previous  iterations,  in 
connection  with  previous  data  items  in  the  series,  replacing  the 
segments  Sj  in  the  respective  hypotheses  with  updated  segments  Sj' 
(step  106) .  In  addition,  the  data  segmentation  module  20 
establishes  two  new  hypotheses,  one  hypothesis  comprised  of 
the  original  segments  and  the  other  hypothesis  Hjn  comprised  of 
the  orignal  segments  plus  a  new  segment  Sjq  representing  the  new 
data  item.  The  hypothesis  Hifa,  since  it  contains  only  the 
original,  non-updated  segments,  represents  the  likelihood  that 
the  new  data  item  is  a  "false  alarm",  that  is,  that  it  neither 

l 

j  belongs  to  any  segment  Sj  nor  is  the  first  data  item  of  a  new 

I  segment  S^,.  The  hypothesis  11;^,  on  the  other  hand,  represents  the 

I 

:  ! 

i|  likelihood  that  the  new  data  item  is  the  first  data  item  of  a  new 

!  ! 

;  j 

I  segment  S^,  and  that  the  other  segments  Sj  are  incorrect 
hypotheses. 

The  data  segmentation  module  20  then  proceeds  to  a  series  of 
steps  to  generate  several  likelihood  scores  for  each  hypothesis. 
In  particular,  the  data  segmentation  module  generates  a  raw 
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li]^0lihood.  score  P(Hj/n+l)  for  the  original  (non  updated) 
hypotheses  H;  as 

1  - 

1  where  "Pn"  represents  the  a  priori  likelihood  that  the  data  item 
I  Zi(t)  starts  a  new  segment,  "Pfa"  represents  the  a  priori 
;■  likelihood  that  data  item  z,(t)  is  a  "false  alarm,"  that  is,  that 

,  j 

i  it  does  not  belong  to  any  segment,  and  "K"  is  the  number  of 
segments  Sj  in  the  collection  of  segments  in  hypothesis  Hj.  The  a 
priori  likelihood  values  are  provided  by  the  a  priori  knowledge 

i 

input  module  21,  and  are  generated  in  any  conventional  manner. 

The  data  segmentation  module  20  then,  for  each  hypothesis  Hj  in 
the  collection  of  hypotheses  H  updated  in  step  106,  generates 
likelihood  scores  for  a  series  of  hypotheses  Hjj,  where  each 
hypothesis  Hy  corresponds  to  the  collection  of  segments  in 

i 

hypothesis  H;,  but  replacing  the  original  of  segment  Sj  with  the 
updated  segment  Sj',  as  well  as  for  the  hypotheses  and  Hjfa 
I  (step  110) .  In  those  operations,  the  data  segmentation  module  20 
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the  likelihood  score  P(H2^)  for  the  augmented  hypothesis  (that 
is,  the  hypothesis  that  the  data  item  is  the  first  data  item  for 
a  new  segment)  as: 

k 

=  P{H./n+l)  n 

I 

I  ’ 

ii  and  the  likelihood  score  P(HiFA)  false-alarm  hypothesis  Hipa  as 

Pi^iPA)  =  PiHi/n  +  l)PpA 

After  generating  the  likelihood  scores,  the  data 
segmentation  module  2  0  prunes  the  hypotheses  Hjj,  HjN  and  HjjpA  ty 
deleting  the  hypotheses  that  have  likelihood  scores  below  a 
predetermined  threshold  value  (step  111) .  The  data  segmentation 
module  prunes  a  segment  Sj,  that  is,  it  completely  eliminates  the 
segment,  when  the  segment  is  no  longer  contained  in  any 
hypothesis  Hy  for  any  index  "i”. 

After  performing  steps  100  through  111  for  one  data  item, 

i 

I I  the  data  segmentation  module  20  returns  to  step  100  to  process 

■  !  I 

il  the  next  data  item.  The  data  segmentation  module  20  performs 

1  steps  100  through  111  for  each  data  item  representing  a  signal  it 
receives  from  the  measurement  track  formation  module  16. 

After  generating  the  hypotheses  Hj  for  a  set  of  data  items, 
the  data  segmentation  module  20  then  transfers  the  set  of  pruned 
segments  Sj  contained  in  the  hypotheses  Hy  and  Hjn  to  the 
trajectory  estimation  module  22.  The  trajectory  estimation 

14 


1 

module  22  then  performs  a  discrete  grid  search  procedure  depicted 

2 

in  FIG.  3  in  connection  with  all  of  the  segments  Sj  to  select  one 

3 

segment  Sj  as  being  most  representative  of  the  information 

4 

j^epresented  by  the  data  items.  With  reference  to  FIG.  3,  the 

5 

trajectory  estimation  module  22  uses  as  the  various  target 

6 

variables  representing  the  target  states  such  target  variables  as 

7 

range  "r",  bearing  "b",  speed  "s",  course  "c”  and  course  rate 

8 

!  j 

"c",  and  establishes  a  series  of  "bins"  with  the  minimum  and 

9  1 

i 

i  maximum  values  for  each  of  these  variables  as  determined  from  a 

1 

10 

priori  knowledge  of  the  possible  target  (step  120) .  The  result, 

11 

1 

if  the  target  variables  r,  b,  s,  c  and  c  are  considered  to  form 

12 

a  five-dimensional  space,  is  a  five-dimensional  grid  of  a  size 

13 

determined  by  the  minimum  and  maximum  values  for  each  variable. 

14 

In  selecting  a  segment  Sj  as  the  potentially  correct 

15 

segment,  the  trajectory  estimation  module  22  makes  the 

16  1 

j 

determination  based  on  certain  ones  of  the  variables,  as 

17 

indicated  by  the  nature  of  the  particular  data  items,  in  this 

18 

lease  range  "r"  and  bearing  "b" .  After  the  grid  is  established. 

i 

19  1 

I f or  each  discrete  point  in  the  grid,  the  trajectory  estimation 

( 

20 

module  22  generates  a  marginal  density  value  Pjj  (step  121)  along 

21 

the  coordinates  "i"  and  "j"  of  these  data  items  as 

22 

p..  =  x;  e  2  (16). 
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The  trajectory  estimation  module  22  then  identifies  the  point 
(i,j)  at  which  the  marginal  density  value  is  a  minimum  (step 
122) ,  and  adjusts  the  perspective  of  the  grid  so  that  it  is 
centered  over  that  point  (123).  The  trajectory  estimation  module 
22  determines  whether  a  selected  accuracy  level,  as  determined  by 
the  resolution  of  the  grid  generated  by  the  trajectory  estimation 
I  module  22,  has  been  reached  (step  124),  and  if  not  returns  to 
li  step  12  0  to  repeat  the  operations.  The  trajectory  estimation 

i 

j 

module  22  repeats  the  operations  until  the  selected  accuracy 
level  has  been  reached,  and  if  so  it  performs  a  chi-squared  test 
in  connection  with  the  segment  Sj.  It  will  be  appreciated  that 
the  trajectory  estimation  module  22  will  perform  these  operations 
in  connection  with  all  of  the  segments  Sj,  and  it  will  use  the 
results  of  the  chi-squared  tests  for  all  of  the  segments  to 

I 

identify  one  as  being  most  representative  of  the  data,  and  that 
segment  is  selected  as  the  correct  one. 

As  noted  above,  and  returning  to  FIG.  IB,  the  segment 
selected  by  the  trajectory  estimation  module  22  is  coupled  to  a 

1 

trajectory  analysis  and  validation  module  23  and  a  trajectory 
I  characteristics  module  24  for  verification  using  conventional 
i  statistical  measures  testing  the  likelihood  or  probability  that  a 
trajectory  is  representative  of  the  information  contained  in  the 
signals  received  by  the  sensor  arrangement  11. 

The  invention  provides  a  number  of  advantages.  It 
facilitates  the  detection  and  use  of  use  of  signals  received  that 
are  below  the  (initial)  signal-to-noise  ratio,  despite  the  fact 
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that  such  signals  are  embedded  in  an  increased  level  of  clutter 
and  noise.  This  enables  the  system  to  determine  the  location  and 
bearing  of  a  target  relatively  early  and  quickly. 

In  addition,  the  invention  provides  an  arrangement  for 
quickly  and  reliably  tracking  a  target  with  an  environmentally 
perturbed  minimal  data  set  comprising  only  a  few  data  items  which 
i|  exploits  target  kinematics  and  a  priori  knowledge,  which  further 

i ! 

n  allows  for  anomolies  in  the  data.  Using  multiple  hypothesis 
i|  techniques,  such  as  as  described  in  connection  with  the  data 

i  ‘ 

I  segmentation  module  20  and  the  trajectory  estimation  module  22 

i 

(FIGS.  2A,  2B  and  3)  allows  for  accommodation  of  changes  in 
dynamics  and  quik  evaluation  of  the  target  dynamics.  The 
segmentation  of  acoustical  information,  as  performed  by  the  data 
segmentaiton  module  20,  allows  for  partitioning  of  the  data 
according  to  similar  features,  which,  in  turn,  allows  for  rapid 
detection  of  motion  changes.  The  discrete  grid  search  technique 
I  performed  by  the  trajectory  estimation  module  22  provides  for 
I  relative  stability  in  the  non-linear  estimation  problem  and 
I  exploitation  of  a  priori  knowledge  of  target  motion, 
i  It  will  be  appreciated  by  those  skilled  in  the  art  that  the 

■I 

i;  new  arrangement  can  be  implemented  using  special-purpose  hardware 

i! 

j  or  a  suitably-programmed  general  purpose  computer. 

I  The  preceding  description  has  been  limited  to  a  specific 

embodiment  of  this  invention.  It  will  be  apparent,  however,  that 
variations  and  modifications  may  be  made  to  the  invention,  with 
the  attainment  of  some  or  all  of  the  advantages  of  the  invention. 
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Therefore,  it  is  the  object 

such  variations  and  modifications  as  come 

and  scope  of  the  invention. 


t^o  cover  all 
within  the  true  spirit 
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SYSTEM  AND  METHOD  FOR  RAPIDLY  TRACKING 

4 

HIGHLY  DYNAMIC  VEHICLES 

5 

6 

ABSTRACT  OF  THE  DISCLOSURE 

7 

A  trajectory  estimation  system  for  estimating  a  trajectory 

i 

8 

of  a  target  in  response  to  a  series  of  data  items  which  generated 

9  i 

in  response  to  motion  of  the  target.  The  trajectory  estimation 

10 

system  includes  a  data  segmentation  means  and  a  trajectory 

11 

selection  means.  The  data  segmentation  means  processes  the 

12 

series  of  data  items  in  accordance  with  a  regression/multiple- 

13 

hypothesis  methodology  to  generate  a  plurality  of  segments,  each 

14 

having  associated  data  items  which  have  similar  features.  The 

15 

trajectory  selection  means  for  processing  said  segments  in 

16 

accordance  with  a  multiple-model  hypothesis  methodology  to 

17 

generate  a  corresponding  statistically-supportable  candidate 

18 

trajectory  motion  estimate  of  target  motion  thereby  to  provide 

19 

indicia  of  an  overall  trajectory  of  the  target. 

KNOWLEDGE  INPUT 


TO  FIG.  2B 


120.  ESTABLISH  A  GRID  USING  PARTITION 
OF  TARGET  STATES  WITH  SPECIFIED 
MINIMA  AND  MAXIMA  USING  AVAILABLE 
KNOWLEDGE  OF  TARGET  BEHAVIOR 


121.  FOR  EACH  POINT  IN  ESTABLISHED 
GRID,  GENERATE  A  MARGINAL  DENSITY 
PROBABILITY  VALUE 


122.  IDENTIFY  GRID  POINT  HAVING 
MINIMUM  ALONG  AXES  OF  GRID 
ASSOCIATED  WITH  PARAMETERS  OF 
INTEREST 


123.  ADJUST  BOUNDS  OF  GRID  CENTERED 
AT  THE  MINIMUM  POINT 


124.  DETERMINE  WHETHER  DESIRED 
ACCURACY  HAS  BEEN  ACHIEVED 


